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The thermodynamics framework of an interacting quantum gas trapped by an arbitrary external 

' r-| ' 

Q ■ potential is reviewed. We show that for each confining potential, in the thermodynamic limit, there 

^ . 

^ , emerge "generalized" volume and pressure variables V and V, that replace the usual volume and 

hydrostatic pressure of a uniform system. This scheme is validated with the derivation of the virial 

■ expansion of the grand potential. We show that this approach yields experimentally amenable 
procedures to find the equation of state of the fluid, V = V(y /N, T) with the number of atoms, 

, as well as its heat capacity at constant generalized volume Cy = C\;{V, N,T). With these two 

Q ' functions, all the thermodynamics properties of the system may be found. As specific examples 

^ ■ 

we study weakly interacting Bose gases trapped by harmonic and by linear quadrupolar potentials 
, within the Hartree-Fock approximation. Comparisons with experimental results of a ^^Na ultracold 

gas are also presented. We claim that this route should provide an additional and useful tool to 
T^lj- ■ analyze both the thermodynamic variables of a trapped gas as well as its elementary excitations. 
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I. INTRODUCTION. 

The realization of quantum ultracold trapped gases 
has opened new test grounds of both known and novel states of matter. This, in turn, 
has vig orously stimulated first principles analysis of the physics of many body systems |12l. 

25l . |26| . |27|. Among the many fundamental 
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aspects of these fiuids, a ubiquitous characteristic stands out, namely, the fact that these 
systems are spatially inhomogeneous due to the confining non-uniform external potential of 
the trap. Hence, usual theoretical tools, concepts, and the concomitant physical intuition 
that has been developed for uniform systems confined by rigid-wall vessels, must be adapted 
or completely reformulated to deal with the ensuing nonuniformities and the presence of the 
trap itself. The main purpose of the present article, in the light of this fact, is to review 
the emergence of "new" and proper mechanical variables for confined fiuids that modify the 
costumary way one is use to deal with their thermodynamic properties. 

The main differences on which we base our discussion, between the thermodynamics of 
non-uniform trapped gases and that of uniform systems, is that while the hydrostatic pres- 
sure p of the fluid becomes a local quantity, the volume V that the gas occupies is no longer 
precisely deflned. That is, these two variables are not anymore the appropriate thermody- 
namic variables to describe the mechanical equilibrium of the inhomogenous trapped fluid. 
The role of the usual p and V, however, is replaced by a unique pair of variables for each 
confining potential, called generalized pressure V and volume V. With the aid of these 
two variables, thermodynamics can then be fully used and exploited as the usual tool to 
characterize and describe the state of the gas. It is our contention that useful and typical 
properties of fiuids, such as equations of state, heat capacities and compressibilities, have 
not been introduced in the description of the physics of the trapped gases so far, because 
of the lack of the appropriate thermodynamic variables. It is our goal to help filling this 
gap. As we shall discuss and show here, these properties can be extracted from the knowl- 
edge of the density profile and other simple properties of the trapped gases. We emphasize 
that thermodynamics, besides characterizing a macroscopic system, also yields information 
regarding both its microscopic interactions and/or its elementary excitations. 

The present article builds in previous works and results dealing with a quantum Bose 
gas trapped in a harmonic potential 1^, 



161, 



13, 



181 ] and in comparisons with experiments in 
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ultracold Sodium gases 28|, |29|, [SOj. Here, we extend the analysis to any arbitrary external 



confining potential and show specific results for Bose gases trapped in both harmonic and 
linear quadrupolar traps. 

We shall proceed in the following way. First, we describe the general system under study 
- an interacting fluid confined by an external, inhomogeneous potential Vexti'f^) - and we 
provide several specific examples to introduce the generalized volume V and pressure V. We 
discuss their physical interpretation and the corresponding thermodynamic limit in which 
these variables emerge. We briefly mention additional thermodynamic variables that may 
arise from the external potential itself and from the possibility of varying the interatomic 
interactions externally; the latter is now a reality being widely exploited We 
show how the equation of state, V = V{V/N,T), and the heat capacity at constant volume 
V, Cv = C\;(y,N,T), can be measured within the current experimental setups. We then 
explictly deal with a confined quantum gas including interatomic interactions, and we derive 
the virial expansion of the grand potential to validite our approach. As a corollary of 
our treatment we show that the so called "local density approximation" (LDA) follows 



within the appropriate thermodynamic limit of the confining potentials [3l|, l32|, |33| . Because 
the density inhomogeneities appear at macroscopic length scales, it may be expected that 
LDA should apply, nevertheless we do not make any explicit assumption on the density 
profiles and it is thus reassuring to find that the virial expansion when applied to confined 
quantum fiuids yields LDA. As we shall discuss, however, care must be taken when using 
it. We shall see that while it gives a procedure to calculate thermodynamic properties of a 
confined inhomogeneous fiuid, it does not imply that the local states of the trapped fiuid 
are thermodynamic states of the corresponding homogenous fiuid. 

We then devout a section to apply our general framework to a confined weakly interacting 
Bose gas. We study a gas confined both by a harmonic and a linear quadrupolar traps 
within the Hartree-Fock approximation. We calculate the phase diagram V — T and the 
heat capacity at constant volume. A brief discussion on the nature of the normal gas to 
superfiuid transition is presented. We show then a comparison of these predictions with 



experimental results of Bagnato et al. on Sodium gases 



29]. 
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II. THERMODYNAMIC VARIABLES OF TRAPPED GASES 



The system consists of identical atoms or particles of mass m with Hamiltonian 

Hn = T.^ + I: + E Kx.(r1). (1) 

i=l -^"^ i<j i=l 

We assume additive pairwise potentials but the analysis may be extended to arbitrary in- 
teratomic interactions. Below, we shall consider the second quantized version of the above 
Hamiltonian to explicitly include the scattering length. The external potential of the trap 
Vexti'^ confines the system. To serve this purpose, it should have at least one minimum 
and must obey that Vexti'f^ — > oo for \ f\ — oo. For rigid- wall containers it is costumary not 
to write down the potential. Here, we include it as Vexti^^ = if r is within the volume 
V enclosed by the rigid walls and Vexti'f^ = cxd if r is outside of it. Typical examples of 
traps of atomic gases are Vextir) = (l/2)m(u; ■ r)^ a harmonic potential, such as in Ref.jli, 
and Vext{r) = \A ■ r\ a linear quadrupo 



potential such Poschl- Teller 
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ar potential 281] ; but one can consider any confining 
35l | Vextif) = K)/ cos(7 ■ r). This last case serves as an 
example of a potential that while it gives rise to a generalized volume, it also introduces 
an additional intensive variable, namely Vq. We write these potentials to exemplify the 
appropriate thermodynamic variables as well as the thermodynamic limit for each case. 

To illustrate how the generalized variables emerge and how the thermodynamic limit is 
to be taken, we shall deal here first with a classical ideal gas. Below, we shall treat an 
interacting fluid and we shall verify the correctness of the identification of the variables 
given here. Consider, therefore, a system given by the Hamiltonian ([1]) with no interatomic 
interactions, i.e. -udrjjl) = 0. Assume the system is in thermodynamic equilibrium at 



temperature T. The canonical partition function is, 
Z(T,iV, V,r^) = j d'^'p J d'^'r exp 



(2) 



where (3 = 1/kT. We assume the external confining potential to be of the form Vext = 
Vexti^/lxiU l^yi z/lzi^) where the quantities U do not necessarily have units of length and rj 
stand for other parameters, such as Vq in the Poschl- Teller potential above. Integration of 
the partition function yields 



Z{T,N,V,n) = j^m.v)V)\ (3) 



where At = /i/(27rmfcT)^/^ is de Broglie thermal wavelength, V = Ixlyh is the generalized 
volume and the function ({P^v) is defined by 

((P,n)V= [ e-^^-*('')rfV (4) 



Helmholtz free energy is found with F = — kT In Z and, after taking the limit — > oo, 
yields 



F{N,T,V, 7]) = -NkTiln 



+ 1 • (5) 



For the free energy per particle, F/N, to remain finite in the thermodynamic limit, — > oo, 
it must be required that the "generalized" volume diverges, i.e. V oo, keeping constant 
the "density" N/V. As it will fully justified below, V is an extensive thermodynamic vari- 
able. The generalized volume certainly is proportional to the actual average volume that 
the system occupies, V ~ ({f3,ri)V. This average volume, however, is not an independent 
thermodynamic variable since it depends on the temperature. Moreover, it is not a correct 
variable since the actual volume that the system occupies is, in general, unbounded. Nev- 
ertheless, the thermodynamic limit V — > oo indeed implies that the volume of the system 
becomes arbitrarily large. 

For the particular external potentials here considered one finds, V = V for rigid walls, 
V = 1/uJxUJyUJz for the harmonic potential, V = l/AxAyA^ for the quadrupolar potential, 
and V = l/7x7j/7z for the Poschl- Teller potential. Likewise, we find = 1 for rigid 
walls, (,{(3) = (27r//?m)^/^ for the harmonic potential, ((P) = Svr//?^ for the quadrupolar 
potential, and ({l3,Vo) = 47r /^^^2 exp[— /9Vo/ cosx]dx for the Poschl- Teller potential. For 
the harmonic case, it has been known for quite a long time that the thermodynamic limit is 



the one here presented 361]. As we shall see, while the role of the generalized volume is com- 
pletely analogous to that of the usual volume in homogeneous systems, the thermodynamic 
properties of the different confined fiuids show very strong variations on their temperature 
dependences due to the function ({P,ri). 

As it will be fully justified below with the virial expansion for an interacting gas confined 
in an arbitrary potential, the generalized volume is a bona-fide extensive variable. Therefore, 
there exists an intensive variable, conjugate to the volume V, that we call the generalized 
pressure V and given by 



N,T 
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Here, F = F{N,T,V) is Helmoltz free energy including interatomic interactions. By a 
simple calculation one obtains, 



1 



exty i) 

1=1 



= ^ / p(r-')f-VKx.t(r)rfV, (7) 

where in the first line the average is performed in the corresponding ensemble and in the 
second line we have introduced the density profile p(r). The last equality is a very useful tool, 
as we shall insist throughout this paper: It gives a direct way to calculate the generalized 
pressure with the sole knowledge of the external potential Vexti'f^, the density profile p{r) 
and the temperature T. That is, the equation of state V = ViV /N,T) is a measurable 
quantity. In Section V we shall use this result for the calculation of the phase diagram of a 
confined interacting Bose gas. 

It is important to point out that the identification of the generalized pressure is not only 
a formal one but it has a clear physical meaning. From the thermodynamic definition of 
V, Eqs.(l6]) and ([7j), one sees that the product VV equals (1/3) of (minus) the virial of the 
external force. Hence, one can recall that mechanical equilibrium in a fiuid is given by Pascal 
law, 

V-P(r) = -p(r)VK.t(r1, (8) 

where P{f) is the pressure tensor of the fiuid. One expects the pressure tensor to be a local 
quantity, P(r) = p{r)l, where 1 is the unit tensor and p{r) the local hydrostatic pressure, 
barring phase-separated states within the confined fiuid. By calculating the virial of the 
right hand side of Eq.Q, and after integrating by parts, one finds 

J TrP d^r = J d^r p{r) f ■ VVext{r) 

= 3VV. (9) 

That is, V for a non-uniform fiuid confined by a given external potential plays the same 
thermodynamic role as the hydrostatic pressure p in a uniform fiuid: it is the quantity that 
bears the information that the fiuid is in mechanical equilibrium. It is somewhat puzzling 
to realize that expression yields only an identity for the rigid-wall case and does not give 
a calculational tool for the hydrostatic pressure. The latter needs the knowledge, at least 



for pairwise interatomic interactions, of the two-body density correlation function 37|]. Here, 



we find that for inhomogeneous systems knowledge of the one-body density suffices. 



We note that with the above identification the change in free energy is 



dF = -SdT - VdV + ^idN 



(10) 



with S the entropy and /i the chemical potential, provided the rest of the intensive parameters 
of the external potential, as well as those of the interatomic potential u{rij), see below, 
remain constant. For instance, if the external potential has additional intensive parameters 
7], such as Vo in the Poschl- Teller potential, and these are also externally modified, there are 
additionally changes in the free energy given by 



35] 



Tdrj, 



where F stands for conjugate extensive variables to t], 



(11) 



(12) 



T,N,V 



The variables rj enter through the function ({jS, rj) only, see Eq.Q. From the above relation- 
ship one finds that an external potential can give rise to as many thermodynamic variables 
(and their conjugates) as the number of parameters needed to specify it. However, there is 
only one variable that plays the role of the volume. For simplicity, we shall assume that the 
additional intensive variables r] of the external potential remain constant. 

The measurement of the heat capacity Cy should also be achievable within the current 
experimental setups. The proposal of this measurement consists of an adiabatic compression 
or expansion. We note first that the ultracold trapped gases are actually isolated and 
confined by magnetic or optical traps [ll, 3, [sl, 4, 0, [a, Q, [sl, [q, 10, 11|. Therefore, a slow 
change of the confining potential, namely, of the generalized volume V, should give rise to 
an increase or decrease of the temperature T depending on whether the fluid is compressed 
or expanded. Since the generalized volume and temperature are measurable in the current 
experiments, the quantity {dV / dT)^ s can, therefore, be calculated. The corresponding heat 
capacity can then be evaluated using then the following thermodynamic relationship. 



dV 



(13) 



We note that previous knowledge of the equation of state is needed for the calculation of the 
second factor on the right hand side of (fT3|) . However, the measurements of the equation of 
state and of the quantity {dV /dT)isi^s correspond to two different sets of experiments. 



It is also of interest here to mention that one of the most interesting and novel aspects 
of the current ultracold gases is the fact that the interatomic interaction potential may be 
modified externally by means of magnetic fields , |ll|]. This, in turn, traduces into 



an external modulation of the scattering length. That is, the scattering length becomes a 
(intensive) thermodynamic variable itself and, therefore, there appears a thermodynamic 
extensive variable conjugated to it. To be explicit, we write down the Hamiltonian in a 
second quantized version and consider a contact interatomic potential. 



H 



(14) 



jklm 



where n and e„ stand for the eigenstates and eigenvalues of the three dimensional one-particle 
Hamiltonian in the presence of the external potential, aj^ and a are creation and annihilation 
operators. The "prime" in the second sum refers to the restrictions introduced by assuming 
an isotropic two-body potential. The coupling parameter is U = A-Kfi^a/m^ with m the 



atom mass and a the scattering length 
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the relevant quantity is rather than U 



38 



:\ has already been recognized, however, that 



39l | . Therefore, since U is clearly an intensive 



quantity, there exists an extensive variable C given by. 



^ " [dm) 



T,N,V,'q 



iklm 



(15) 



This quantity has been called the contactp^. |40| and it is also related to the large momentum 
tail of the particles momentum distribution. Thus, this yields an additional term in the 
change of free energy given by 



(16) 



In the current literature C has been introduced as the adiabatic change of the internal energy 
E with respect to l/f/[39^. From the previous analysis, one finds that such a definition is 
actually incorrect. Since C is the isothermal change of the free energy with respect to 1/U, 
standard thermodynamics 4l|] indicates that in the microcanonical ensemble the correct 
relationship is 

1 /f)F.\ 

(17) 



1 

U 



'dC 



s,N,v,r 



A full study of the variables C and 1/U is of fundamental importance, but its is out of the 
scope of the present article. Nevertheless, as an instance of the relevance of the mechanical 



8 



variables with these studies, we mention here the recent interest on the behavior of Fermi 
gases (^°K and ^Li) near the unitarity hmit where the scattering length diverges, i.e. the limit 
1/U — 0, because it appears that thermodynamics becomes universal, that is, independent 
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25 



26 



271]. In this region, the 



of the interatomic interactions; see e.g. Refs.0,y,y,y,i 

gas behaves as an ideal one in the sense that it obeys the ideal virial relationship, which for 
a harmonic trap reads E = 2N < Vext >, with E the internal energy. The connection to 



the present work is that < Vext > 
as ideal or not, see Eq.© and Ref. 



Refs 



or a harmonic trap is 3/2VV, whether the gas behaves 



26 



ISj. That is, the quantity that has been measured in 



27| is precisely the generalized pressure for the harmonic trap. And indeed, if 



the universality hipothesis is correct 25|| , using the virial theorem for ideal gases for arbitrary 
potentials, yields the following relation that should be obeyed in the unitarity region. 



N o 

E=<Y.Vext{r^>+-VV. 
1=1 ^ 



(18) 



Away from the unitarity limit, this equation is no longer valid but the measurement of P, and 
of the heat capacity as we described above, can be performed to obtain the thermodynamics 
of those states. 



III. VIRIAL EXPANSION FOR ARBITRARY CONFINING POTENTIALS. 

With the identification of the generalized variables and the corresponding thermodynamic 
limit in hand, we now turn to the virial expansion of an interacting quantum gas. The 
classical case appears as the limit of high temperatures, as we shall indicate it. We extend 



the analysis described in advanced textbooks on statistical physics |42l. I43l . I44| . With this 
expansion we shall validate the thermodynamic variables introduced in the previous section. 
Again, we assume the system is in thermodynamic equilibrium at temperature T and we 
analyze it in the grand canonical ensemble. We thus consider a chemical potential yU whose 
value may be found by imposing a given number of particles A^. The grand potential is given 
by, 

oo 

= -kT In e'''''^ Tr'e-^^^, (19) 
Ar=0 

where 

Tr' e-/^^- = ^ 1 d^'^T W^^(fi, r2, . . . , r^) (20) 
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and 

Wwin, fa, . . . , rW) = ^ < n, fa, ... , f7v|e"^^^|fip, fap, . . . , f^vp > . (21) 

p 

The sum is over all permutations of 1, 2, and e = ±1 for bosons or fermions. 



To find the virial expansion, equation fll9p is first rewritten as, 421. l43l. |44 1 



OO 1 

-pn = J2 e^^" -, In, (22) 



n=l 



where the functions J„ are given by 

In = J d?''rUn{n,...,rn), (23) 
and, in turn, the Ursell functions are found from the hierarchy, first order, 

f/i(l) = 1^1(1), (24) 

second order, 

f/2(l,2) = 1^2(1,2) - f/i(l)t/i(2), (25) 

third order, 

t/3(l,2,3) = l^3(l,2,3)-t/i(l)f/2(2,3)-f/i(2)f/2(l,3)-f/i(3)f/2(l,2)- 

f/i(l)[/i(2)f/i(3). (26) 

and so on. 

The virial expansion follows after finding the value of each contribution J„ in the thermo- 
dynamic limit. We do this now for a general confining external potential Kxi(^)- We proceed 
by systematically calculating /„ order by order and then generalize it to /„. We have done 
so from Ii to 1^. Although lengthy, the corresponding calculations are straightforward and 
we explicitly present in the Appendix the case I2 only. Next we discuss the results. 

The calculation of Ji is very simple but serves to indicate an important point: 

Ii = j d^r < f\e~^"'\f>, (27) 
where the one-particle Hamiltonian Hi is given by 

Hi = ^ + Ve,t{f^. (28) 
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In the thermodynamic hmit, equivalent to the use of the so-called semi classical density of 
states[4^, one readily finds 

h = ^C(/9)V. (29) 

The same expression would be obtained if the system obeyed classical mechanics. As we 
indicate below, this result is always obtained for the center of mass motion of the n-particle 
cluster involved in /„. 

Following the explicit calculation of I2 shown in the Appendix suffices to see how to find 
In- The key is in the separation of center of mass and relative coordinates. This change 
of variables is generally, R = ^{fi + r 2 + . . . + fn), f^^^ = ^1 — ^2? f^'^^ = ''^2 ~ "^s? • • 
r^"~^^ = r„_i — fn, with their canonical conjugate momenta. The main assumption is that, 
in the thermodynamic limit V 00, the following approximation is correct, 

Kxi(ri) + Kxt(r2) + . . . + Kx.t(r„) ~ nVe.t{R), (30) 

where (ri, r2, . . . , fn) are to be given in terms of the variables {R, f^^\f^'^\ . . . , f("^i)) by the 
above transformation. This approximation separates the center of mass motion from the 
relative ones of the corresponding n-particle cluster. The ensuing motion of the center of 
mass is always quasiclassical and its contribution to /„ is proportional to ({nP)V / X^, see 
Eq.f l29l) . Within the same limit, the contribution from the relative coordinates yields the 
usual quantum virial coefficients bn{T), the same as those calculated in the uniform case, 
i.e. they are universal for all confining potentials. Thus, one generally finds 

V 



Xrp 



,J{n(3)bn{T). (31) 



The validity of the above procedure requires, first of all, that the interatomic interactions 
must be "short-range" , namely, decaying faster than otherwise the virial coefficients 



b„ do not exist 



31 



41 



42j. For high temperatures, in the classical regime, one finds that the 
intermolecular potential must vanish for lengths r ^ cr, with a the range of such a potential. 
At low temperatures the bound is set up by either the thermal de Broglie wavelength or the 
scattering length a 46|. If the gas behaves as an ideal one, the relevant length is de Broglie 
wavelength. In any case, as long as the relative coordinates remain bounded by a finite 
quantity, however large, one can take the limit of very large volumes V 00 and separate 
the motion of the center of mass from the relative motions of the involved n particles. It 
is of pedagogical interest to mention that the present approximations are always tacitly 

11 



performed in the uniform case of rigid-wall potentials, by neglecting boundary terms after 
the corresponding change of variables. 

Summarizing, we find that in the thermodynamic limit the grand potential can be written 
in general as, 

^=-kTVY.-r^-^UT) (32) 

n=l 

with 61 = 1. This expression for the grand potential validates the introduction of the 
mechanical variables V and V, as we now verify. The number of particles and the entropy 
S can be calculated from (minus) the partial derivatives of Q with respect to fi and T 
respectively. Q, N and S are found to be homogeneous first order functions of V, and this 
implies that V must be an extensive variable and justifies the thermodynamic limit as used 
above. Since it should be obeyed that the conjugate variable isV = —{dVLj dV)T,^, it follows 
that Vt = —VV, as it should. Thus, the generalized pressure is read off (1321) . 

It is interesting to note that the most important difference of the grand potential between 
a given arbitrary external potential and the homogeneous case is the function C{l3) rather 
than the generalized volume V. The latter enters in the same way for any potential, including 
the rigid-walls case; that is, it gives rise to the intensive quantities formed between the 
extensive variables A^, S, E, etc. and V, that remain finite in the thermodynamic limit, i.e. 
N/V, S/V, E/V, etc. However, as it is now well established, the temperature dependence of 
the thermodynamic variables is very different and unique for each external potential. The 
function ({P) gives rise to those differences. 

To illustrate the use of Eq. fl32l) . we apply it to an ideal quantum gas. From the analysis 
in the Appendix and their corresponding value for third and fourth orders, one finds that 
the quantum ideal virial coefficients are given by 

"r-^-^'^^^'"-". (33) 

Hence, the grand potential for an ideal quantum gas can be written as 

-/5^=3TEe^^"C(n/3)-^. (34) 

"^T n=l ^ 

This formula can be directly compared with the corresponding ones for, say, the rigid 
walls potential V = V and Cij^l^) = 1, or the harmonic potential V = l/u^ and 



C(n/3) = (27rA;T/nm)^/^. The "textbook" formulae for these potentials are. [47 1 



Ajf(572)^ e-/^^ - e 
12 



'kT\^ 1 r'^ x-^dx 



for rigid walls, and 

for a 3D isotropic harmonic potential. Expansion of the integrals of these last two equations 
in powers of e^^ yield the virial expansion, Eg. (1341) . 

For completeness of our presentation, we write down the first few terms of the so-called 
virial expansion of the equation of state V = V{N/V,T) for low densities N/V. This can 
be done by finding = N{fi,T,V) from Eq. flS^ and inverting it term by term to yield 
/i = ^{N/V,T), then, substituting the result into V = V{^,T): 



N N 



2^(0) " 'V \('(j3) ' 3(Hj3) " 7 VV 



(37) 

Once again, we remark that the functions ({(3) make all the difference. Since in some 



instances one can refer the calculation of the virial coefficients to a diagramatic expansion 



42|, 



one finds that the diagrams sum up differently for different potentials. We also recall that 



this type of virial expansion was used in Ref. 
Sodium atoms in a quadrupolar potential. 



28l | to fit experimental data from a gas of 



IV. A NOTE ON THE "LOCAL DENSITY APPROXIMATION". 

A corollary from the virial expansion above, Eq.( l32i) . is that the validity of the "local 
density approximation" follows right away from the corresponding expressions for Q, N and 
S. That is, suppose a uniform gas confined by the rigid-wall external potential, V = V 
and C{ni3) = 1. Define the grand potential per unit volume uj{fi,T) = ^l/V, the number 
of particles per unit volume (particle density) p{fi,T) = N/V and the entropy per unit 
volume s(/i, T) = S/V. Now consider the same system but trapped by an external potential 
Vext{r)- Its thermodynamic properties may then be found by implementing the "local density 
approximation": take u, p and s of the homogenous case and make those functions per 
unit volume to be their "local" densities uj{r), p{r) and s{r) in the presence of the given 
external potential, by replacing the chemical potential p by the "local" chemical potential 
Piocaiif) = A* ~ Vexti/f^- It turns out that integration of uj{f), p{r) and s{r) over all space 
yield the exact expansions for fi, and 5, in the presence of Vext{r), as given by Eq.(|32ll 
and its derivatives. That is, one finds that LDA procedure gives rise to exact results. We 
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recall that the validity of LDA for classical and quantum systems in this limit was rigorously 



proved in Refs. 



31 



32| and [33|, respectively. 



The above description does show that in the thermodynamic limit the system is locally 
homogenous and that "locally" actually means in length scales large compared with those 
of interatomic interactions. It is in this latter connection that LDA is largely used without 
the need of further justification. There is a warning, however, that must be raised when 
using LDA. It may appear that if one is able to find any thermodynamic variable q for a 
homogenous system and express it in terms of the chemical potential fi and temperature 
T, namely q = q{fi,T), its local counterpart when in the presence of an external potential 
Vext{r) is simply g(r) = q{fiiocai{^,T). This, in general, is incorrect; its is strictly justified 
for ^l/V, N/V and S/V only . It is incorrect, for instance, for the internal energy and other 
free energies, except Q, as well as for other functions such as heat capacities. This statement 
can be verified by using expression (132!) for the Grand Potential of a confined fluid. Thus, 
the fact that the system is locally homogeneous does not imply that the local states of 
the confined fluid are thermodynamic states of the corresponding homogeneous system, i.e. 
qiocai{r) ^ q{fiiocai{r),T) in general. 

To verify the previous statement, let us consider the internal energy density e = E/V of a 
homogenous system. From the knowledge of p{= —uj), s and p, the energy may be calculated 
as e = Ts — p + fip. Thus, one can find e = e{fi,T). Then, consider a confining potential 
Vext{r) and implement LDA. Presumably, the internal energy of the trapped fluid would be 
E = J e{fiiocai{'f^) , T) d?r. This, however, is not the correct value of the internal energy of the 
fluid in the trap. The use of Eq.(l32l) allows us to verify it by explicitly calculating E for a 
trapped fluid. That is, an alternative form is to use the Grand Potential for the fluid in 
the trap, as given by Eq. fl32l) . with the appropriate variables V and C(/5)- Then, calculate 
the internal energy of the trap as = TS — VV + pN . One finds a different result. This 
can be checked very simply by considering an ideal classical gas, i.e. use hi = 1 and 6„ = 
for > 1 in Eq. (!32l) . The (incorrect) result by the first route is E = 3NkT/2 for all traps, 
while the (correct) result from the second route is 

It is not difficult to show, using Eq. (l32!) . that the local internal energy density is actually 
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given by 

eiocai{r) = e{fiiocai{r),T) + Vext{r)p{^.iocai{r),T). (39) 

That is, integration of eiocaiif^ yields the correct internal energy of the trapped fluid, while 
integration of ei^niocaiif) )T) does not. This result can be physically rationalized by arguing 
that e(/i, T) considers kinetic and interatomic interactions only and, therefore, the addition 
of the external energy density is required. For other functions the discrepancy is more 
difficult to assess. For instance, one can find expressions for the specific heat at constant 
volume using the homogeneous thermodynamics, say Cy /V = c\/(/i, T). Again, even if eiocai 
is used, LDA expression for cv{lJi'iocaiif)^T) does not integrate to the correct value found 
directly from the grand potential for the non-uniform fiuid T, V), Eg. (1521) . Actually, it 
is not even clear if a local specific heat is a meaningful quantity. That is, the heat capacity 
at constant volume translates locally into an specific heat at constant density; however, 
if a confined fiuid is heated up, keeping the same number of atoms without changing the 
trap, its density profile changes. Hence, it is not possible to keep the local density constant 
while heating the system up. For the overall confined fiuid, heat capacities are certainly 
meaningful but these must be calculated or measured either at constant V or V. 

V. EQUATION OF STATE AND HEAT CAPACITY OF WEAKLY INTERACT- 
ING BOSE GASES. 

The chief importance of correctly identifying the generalized pressure and volume, resides 
in its use as a tool to characterize a given system. For the simple case of a one-component 
gas and for a fixed interatomic interaction, there are only two independent thermodynamic 
variables, say, the temperature T and the molar, or per particle, generalized volume v = 
V/N. Therefore, at least two further functions of these variables should be measured or 
calculated independently, in order to obtain the full thermodynamics of the system. We 
choose the equation of state V = V{v,T) and the specific heat at constant generalized 
volume C\)/N = Cy{v,T). As we have explained in Section II, these two quantities should 
be very easily measured in the current experiments of ultracold gases. 

We now turn our attention to the calculation of the equations of state and heat capac- 
ities of a weakly interacting Bose gas confined in a harmonic and in a linear quadrupolar 
potential, Vextir) = (l/2)m(a; ■ r)^ and Vext{r) = \A-f\ respectively, within the Hartree-Fock 
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approximation 13|, |47|, |48| . This is a self-consistent calculation for the density profiles of the 
thermal and the condensate densities, pth{r) and po{r), that leads to the following set of 



equations 



and 



48| 



Pth[r) 



A ji 



gsdP [p - V,,t{f) - 2Uptf,{f) - 2f/po(r)]) 



po{r) = — (/i - Vext{r) - 2Upth{r}) , 
with the constraint that the number of particles is a given value A^, 



N 



pth{r) d?r + po{r) d^r. 



(40) 



(41) 



(42) 



In Eq.fHOl). (73/2(0;) is the usual Bose function 5'„(a) for n = 3/2. Equation fHTj) is to be 
undertstood valid for values when the right-hand-side is positive or zero. As a matter of 
fact, this is how the normal to Bose-Einstein condensation (BEC) or superfluid transition is 
identified, i.e. given the temperature T, the transition occurs for the value of the chemical 
potential below which the condensate density poir) is different from zero. Equation ( HTi) is 
the Gross-Pitaevskii equation in the thermodynamic limit where the kinetic energy term may 
be safely neglected. The above set of equations suffers essentially from the fact that it does 



2l|. However, it 



not consider the expected Bogoliubov excitations at very low temperatures 
should be fine for temperatures near the transition 48 1. 

As it is clear from the above set equations, their solution yields the density profile p(f^ and 
the chemical potential p for given values of the temperature T and the generalized volume 
V = V/N. The value of generalized pressure V{v, T) is found from Eq.Q and, together with 
p{v, T), one can further find the molar Helmholtz potential / = F/N as f{v, T) = —Vv + p. 
The molar entropy s = {df/dT)^ follows and, therefore, the specific heat = T{ds/dT)y. 
Our results are summarized in Figs. 1 to 4. 

Figures 1 and 2 show a few isochores {v = const) of the equation of state for both external 
potentials, comparing the ideal case with the interacting Hartree-Fock approximation. Note 
that although the generalized pressures are quantitatively different, even with different units, 
their qualitative behavior is essentially the same. We make the following comments. First, 
in the ideal case the transition BEC line indicates that the pressure vanishes as T — 0. 
That is, just as in the uniform case 411], the condensate exerts no pressure. This is clearly 
changed once interactions are included: the pressure of the condensate is no longer zero. 
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and even at T = the interactions give rise to a remnant pressure. Second, the transition 
temperature, different for each isochore, is shifted down in the interacting case with respect 
to the ideal one. This downshift or the transition temperature is in agreement with results 



Figure 2 A 




0.00' ' 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 

T 



17 



Figure 2 B 




of more general theories of trapped Bose gases 



19|. Moreover, this is an effect due solely to 



the interactions and not related to finite size effects [4^. And, lastly, the transition line 
in the interacting shall be further described below, marks a smooth normal gas to 

superfluid transition, different to BEC where discontinuities in the second derivatives of the 
free energy are encountered; in the interacting case, up to second derivatives - and appears 
that to higher order as well - the free energy is continuous. 

In Figures 3A and 3B we show the specific heat as a function of temperature, for a 
given isochore, for both potentials. Again the qualitative behavior is the same. Once more, 
we see that the transition temperature in the interacting case, marked with an arrow, is 
lower than that of the ideal case. But more interestingly, one finds that at the transition 
temperature the specific heat does not show its maximum value, but rather its minimum, 
and that the transition is continuous. Thus, it shows sign of being neither a first order 
nor a critical transition. We understand that the present is a mean-field calculation and, 
as mentioned above, perhaps not the best description of a superfluid; however, mean-field 
theories typically yield incorrect quantitative results but do not change the order of the 
transition. The origin of the continuity of all the thermodynamic properties may be traced 
back to the behavior of the condensate fraction. This is exemplified in Figure 4 where we 
compare the condensate fraction Nq/N of the ideal with the interacting case. Below T^, the 
ideal condensate fraction is Nq/N = 1 — (T/Tc)^^'^ for the uniform case, Nq/N = 1 — (T/Tc)^ 
for the harmonic trap, and Nq/N = 1 — (T/T^)^/^ for the linear quadrupolar potential. Above 
Tc, Nq/N = 0. Thus, the transition in the ideal case has a discontinuity in the derivative. 
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Figure 3 A 
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However, for the interacting shown in Figure 4, there appears that this transition 

is completely smooth, with no discontinuity or singularity at all. Although not shown here, 
the corresponding isothermal compressibility of the trapped gases, = — 1/V(9V/9'P)t, is 
also continuous showing no sign of any critical fluctuations. This is also in agreement with 
general results concerning trapped fluids [l^. 

Although no comprehensive experimental comparisons have been made using the present 
mechanical variables, Bagnato et al. have performed initial tests with a gas of ^^Na atoms 
in a harmonic trap [71, l29|. Table I summarizes the comparison. In those experiments, the 
measured quantities are the parameters of the trap, the particle density, the temperature, 
the number of particles and the scattering length. The experimental generalized pressure 
is calculated using Eq.Q and the theoretical one is obtained with the HF approximation 
described above. Clearly, for temperatures above BEC the theory shows excellent agreement 
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Figure 4 




0.0^ 
0.0 



with experiments. For temperatures below - the last two entries of the Table- the Hartree- 
Fock approximation does not fare very well, as expected. It appears that the role of the 
interactions is still underestimated. An improvement using better theories, that include 

n 

the proper role and statistical description of Bogoliubov modesf2l!] and/or including the 
Popov approximation 13l. l47l|. is certainly desirable. Nevertheless, these initial experiments 
do show the usefulness of the knowledge of the equation of state, not only for purposes of 
characterization, but clearly as an additional tool to learn about the elementary excitations 
of the superfiuid state of the ultracold gases. 



VI. FINAL REMARKS 



The main purpose of this brief review if to emphasize that the use of the generalized 
thermodynamic variables here presented should lead to useful results for the analysis of 
the physics of ultracold trapped gases, by simply following the rules of thermodynamics. 
We have shown that, perhaps with an elaborate and lengthy experimental procedure, both 
the equation of state, V = ViV /N,T), and the heat capacity, Cy = Cv(V/A^, T), can be 
readily measured with the current experimental setups. Knowledge of these two quantities 
suffices to know all the thermodynamics of the corresponding system. As we mention in 
Section II, extensions to other possible thermodynamic variables, such as the externally 
varied scattering length, may be readily incorporated. It is clear to us that for this framework 
to be really useful one needs, first, to change the usual intuition on volume and hydrostatic 
pressure to their generalized counterparts and, second, to provide examples where these 
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WV/mN (exp) 


WV/mN (theo) 




35 


388 


377.9 


0.79 


12 


130 


129.5 


0.53 


11.2 


110 


120.9 


0.41 


4.5 


45.4 


48.6 


0.34 


4.3 


42.4 


46.6 


0.30 


0.9 


22.4 


9.7 


0.16 


0.8 


6.8 


8.6 


0.09 


0.33 


4.2 


3.6 


0.10 


0.07 


0.9 


0.52 


0.076 


0.04 


0.7 


0.09 


0.065 



TABLE I: Comparison of experimental generalized pressure V with a Hartree-Fock (HF) calculation 
for an ultracold Sodium gas, see Eq.([7]). V = 1/uJxOJyUJz with oj^ = 27r x 36.7 Hz, = 27r x 120.8 
Hz, and a;^ = 27r x 159.6 Hzj?]. The mass is m = 23 a. u. for Sodium and the scatterin g le ngth 
is a = 65ao (Bohr radii). is the number of particles, which for all these experiments [29|] was 
from N ~ 10^ to A ~ 10^. The transition BEC temperatures have been calculated within the HF 
approach. 



variables lead to new insights. The latter is mainly a job for theory and in this review we have 
presented the study of a weakly interacting Bose gas within the Hartree-Fock approximation. 
The results are promising and compare well with experiments, but certainly improvement 
is needed specially for very low temperatures. 

We have also shown that the present approach leads also to validate the use of LDA, 
although we have indicated that care must be taken when using it. This should be more 
delicate when dealing with phase-separated fluids where it is not clear if LDA suffices for 
their description since the interfacial widths of the phase boundaries are expected to be of 
the order of the range of the intermolecular interactions 37]. This comment may apply to 



the states found in Fermi trapped gases of ''Li atoms[8|, |9|, [ll], where there is evidence that 
the confined fiuid phase-separates into a superfiuid and a normal paramagnetic gas, showing 
an interfacial phase boundary. In general, for such inhomogeneous states, one should not 
expect a local picture to be valid across the interface; the thermodynamic potentials are 
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expected to be non-local on the density profiles. An approach based on the generalized 
variables here analyzed may lead to a global and novel picture of those states. 
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APPENDIX: DERIVATION OF h 

In this Appendix we provide the derivation of second order contribution I2 to the virial 
expansion of the grand potential, Eg. (1221) . 

/^From the general expressions in Section III, Eqs. (IT^ - (|23|) . one finds, 

h = j d^ri j d^r2 (< fi, f2|e"^^^|fi, f2 > +e < fi, f2|e"^^^|f2, ri >) - 

j d^ri < fi\e~^"^\ri > J rfVa < f2\e'^"'\r2 > (A.l) 

where e = ±1 for bosons or fermions. The Hamiltonians H2 and Hi are given by Eqs.([T]) 
and fl28p . We make the change of variables to center of mass and relative coordinates, 

R={ri + f2)/2 and f=fi — 'r2, (A. 2) 

with their canonical momenta P and p, respectively. This gives, 

H2 = + VUR + ^) + Ve.tiR - ^) + H2, (A.3) 

where 

n2 = + uir) (A.4) 

2[m/2j 

is the two-particle relative coordinate Hamiltonian. The thermodynamic limit consists in 
approximating 

— * — * 

VeAR+l) + ^-^^^^ ~ 2V;.t(^). (A.5) 

since the range of the center of mass motion becomes arbitrarily large as V — > 00, while the 
ensuing relative motion is bounded. It is easy to verify that the latter is bounded by the 
interatomic interaction range a for large temperatures, as the motion is classical, while it is 
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bounded by the thermal de Broghe wavelength or the scattering length, for low temperatures 
when the motion is quantum. 

The above limit, Eq. (lA.5p . separates the motion of the center of mass from that of the 
relative coordinates, 

j dh (< r\e-^^^ |f > +e < rle"^^^ | - f > - < r\e~'^'p' '^^"''^'^\f >) . (A.6) 

We note that the center of mass motion corresponds to a one-particle system of mass nm 
moving in an external potential nVext{R)- In the thermodynamic limit one is allowed to use 
the semiclassical density of states for the center of mass motion, yielding, 

h = ^CmhiT), (A.7) 

where 

VC(2/?) = J d^R e-2/^^-*(^), (A.8) 
and where the quantum second virial coefficient is given by, 

62(r) = 2'/^\l I dh (< rle-^^^lf > + 

e < r\e-^'^^ \-r>-< f^^e'^^''/'^'^/^) |^ >A _ (^.9) 



As a rule, in the thermodynamic limit the center of mass motion is always quasiclassical 41]. 



3e seen to be the correct 



4l|. For slow collisions. 



The expression for the quantum second virial coefficient above can 
one by comparing, for instance, with the expression given in Ref. 
the relevant ones for ultracold gases, 62 depends on the scattering length a and this may 
become quite large near a Feschbach or potential resonance. The formulae here derived may 
then not be applicable very near such a point, called the unitarity limit, but as it has been 



shown 251] this may be expected since in such a limit the system behaves as if near a critical 
point. We add that the description of the scattering al low energies near resonances is valid 
for interatomic potentials u{r) that decay at least as l/r^[46j]. In the classical limit one also 
finds that the interaction must be "short-range", otherwise 62 does not exist. 

To end this part, we find illustrative to calculate 62 for an ideal quantum gas, i.e. for 
u{r) =0. One finds the so-called "exchange" contribution to the second virial coefficient: 

bf=e^,X%. (A.IO) 
23 
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FIGURE CAPTIONS 

Figure 1. Phase diagram P — T of a weakly interacting Bose gas confined by a harmonic 
trap (A) and by a hnear quadrupolar potential (B). Several isochores {V /N = constant) are 
shown in dotted lines. The solid line shows the normal gas to superfluid (BEC) transition. 
Compare with Figure 2, the ideal case. In the latter the BEC transition line occurs at higher 
temperatures than in the interacting case. Note also that in the interacting case, below the 
transition temperature, the condensed phase exerts pressure. Units are h — 1, m — 1 and 
ao = 1. 

Figure 2. Phase diagram P — T of an ideal Bose gas confined by a harmonic trap (A) 
and by a linear quadrupolar potential (B). Several isochores {V/N — constant) are shown 
in dotted lines. The solid line shows the BEC transition. See caption of Figure 1 for further 
details. 

Figure 3 Specific heat at constant generahzed (molar) volume vs. temperature T for 
a gas confined by a a harmonic trap (A) and by a linear quadrupolar potential (B). The 
dotted line is the ideal case and the solid line the weakly interacting Bose gas. The transition 
temperature in the interacting case is marked with an arrow. Note that while in the ideal 
case the heat capacity is discontinuous at the transition temperature, it appears that in the 
interacting case it is completely continuous. Further, in the latter case, the transition does 
not occur at the maximum value of Cy but at its minimum. Units are h — 1, m — 1 and 
ao = 1. 

Figure 4. Condensate fraction Nq/N as a function of temperature T for an interacting 
Bose gas confined by a linear quadrupolar trap. Note that the derivative of the curve is 
continuous at the transition. See text for details. Units are ^ = 1, m = 1 and Qq — 1. 
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